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Thomas D. CoheiE and Ehzabeth S. WerbofQ 
Department of Physics, University of Maryland, College Park, MD 20742-4111 

The response of the QCD vacuum to very large static external magnetic fields {qB ^ -^qcd) 
is studied. In this regime, the magnetization of the QCD vacuum is naturally described via per- 
turbative QCD. Combining pQCD and the Schwinger proper time formalism, we calculate the 
magnetization of the QCD vacuum due to a strong magnetic field at leading order (one-loop) to be 
proportional to BlogB. We show that the leading perturbative correction (two-loop) vanishes. 



I. INTRODUCTION 

Strongly interacting matter is often studied via elec- 
tromagnetic probes. For example, the electromagnetic 
form factors of hadrons give us information about the dis- 
tribution and motion of quarks in the hadronjj]. Other 
useful probes of hadronic — and nuclear — systems include 
real and virtual compton scattering and various inelas- 
tic processes including those in the deep inelastic regime. 
However, it is less commonly appreciated that-at least in 
principle — electromagnetism can probe interesting prop- 
erties of the QCD vacuum. In essence an external electric 
or magnetic field will polarize the QCD vacuum and the 
response of the vacuum to an external quasi-static field 
gives significant information about the QCD vacuum. 

In practice, the effects of an electromagnetic field on 
observables associated with the QCD vacuum are very 
small for fields achievable in the laboratory or for known 
astrophysical sources. For example, the proposed Ex- 
treme Light Infrastructure (ELI) project [3|| would pro- 
duce fields which are several orders of magnitude smaller 
than needed to have significant effects on the QCD vac- 
uum. Similarly, the most intense known extended fields 
in astrophysics — the surface of magnetars — while large 
enough to cause interesting nuclear physics effects (such 
as causing the binding of diproton states Q) are about 
two orders of magnitude too small to be of relevance 
to the QCD vacuum. Thus the question of how the 
QCD vacuum responds to extended quasi-static electric 
or magnetic fields remains largely of theoretical interest. 
Nevertheless the question is of significant interest and 
has been the subject of study nearly two decadesQ. 

Much of the study of the response of the QCD vac- 
uum to external fields has been in the context of "QCD- 
inspired" models rather than QCD itself. These include 
the Nambu and Jona-Lasino (NJL) model[^, 0, 0, Q, the 
linear a model0, [IH, and a response due to meson 
loops with the mesons described in a quark model p^. 
All such approaches suffer from the fact that it is diffi- 
cult to know how well they actually reflect QCD. Apart 
from this general concern, there are deep reasons to sus- 
pect that models such as the NJL model and the linear 



sigma models are unlikely to correctly reproduce QCD 
in interesting regimes. Calculations in these models are 
based on mean-field theory which is only known to be 
justified in the large Nc limit of QCD. However for weak 
fields (and small pion masses) the system is near the chi- 
ral limit; it is well known that the chiral limit and the 
large Nc limits do not in general commute and thus 
it would not be surprising if the models fail to reproduce 
the behavior of QCD for weak fields. As it happens, these 
models fail to reproduce the qualitative feature of QCD 
for weak fields |14|. 

Fortunately, in the weak-field regime where the mean- 
field models are particularly suspect, there is a reliable 
model-independent approach to the problem — namely, 
chiral perturbation theory (^FT) 15]. This approach is 
based on a scale separation between the masses of the 
pseudo-Goldstone bosons and the characteristic hadronic 
scale (^ 1 GeV). The approach is in essence an ex- 
pansion in the ratio of light scales (psuedo-Goldstone 
masses, momenta) to the typical hadronic scale. The 
effective theory contains a fixed number of constants at 
each order in this expansion. In principle, these con- 
stants are all determined from QCD. However, the ap- 
proach has predictive power even when QCD cannot be 
fully solved. By fitting these constants from some finite 
set of experimental observables, one can predict all other 
observables — albeit only approximately with an accuracy 
fixed by the order in the expansion at which one is work- 
ing. In the pre sent context, a crucial observation was 
made in ref. [14| : an external magnetic field is a source 
of chiral symmetry breaking and if the magnitude is small 
(i.e., eH < A^^^^.^^^;^ where Anadronic ~ IGeV is a typi- 
cal hadronic scale or several times Aqcd), its effects can 
be computed in chiral perturbation theory. This work 
considered the limit of massless pions, worked at low- 
est non-trivial in xPT and was restricted to purely mag- 
netic fields. In this regime, the computation is relatively 
straightforward — it amounts to a one-pion loop calcula- 
tion for charged pions with the interactions with the ex- 
ternal field included to all orders. Technically this is im- 
plemented via the Schwinger proper time formalism [l6j. 
Subsequent work extended the analysis to finite mass pi- 
ons and to the situation of electric as well as magnetic 
fields [13 and to next-to- leading order in yPTflsj. 

The small field regime was tractable precisely because 
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chiral perturbation theory breaks down. In this 
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regime one must either directly solve QCD — for example, 
via lattice simulations — or resort to models. The purpose 
of the present paper is to explore the strong field region: 
eH ^ In such a regime, there is again a scale 

separation and one might hope that the this regime is also 
tractable due to a model-independent and systematic ex- 
pansion. Before proceeding it is useful to recall that for 
the foreseeable future the prospects for obtaining fields 
in either laboratory or astrophysical settings which are 
large enough to have noticeable effects on the QCD vac- 
uum remain dim. The prospects of getting to the strong 
field regime where the fields are large on hadronic scales 
is even more remote. 

Thispaper is not the first to focus on this regime: 
ref. [12| considered strong magnetic fields in the con- 
text of a model. The model has two basic assumptions: 

i) that the dominant effects can be extracted from two- 
body interactions between quarks and anti-quarks, and 

ii) these two-body interactions can be well-approximated 
by a potential model. The dynamics of qq interactions are 
altered due to the presence of the magnetic field — which 
forces the quarks into relativistic Landau orbits and re- 
duces the strong interactions into something which is ef- 
fectively one dimensional. At a sufficiently strong field 
this induces condensation of spin polarized pairs, which 
yields a magnetization of the vacuum. In the case of very 
strong fields it is argued that the system will be tightly 
bound — in which case the potential can be well approxi- 
mated by a color coulomb force with a coupling running 
via perturbative QCD. This interaction can be solved in 
the WKB approximation. A rough estimate for the den- 
sity of induced pairs is given: it is simply the density 
at which the (color-singlet) pairs begin to overlap and 
is thus determined by the size of the pairs. With these 
assumptions the magnetization associated with a given 
flavor of quark is given by 
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where m is the constituent quark mass and Nc and Nf 
are the numbers of colors and flavors respectively in the 
theory. 

It is not immediately clear whether t he p redictions of 
ref. are robust. It was noted in ref. [l^ that quanti- 
tatively, the calculation will have non-trivial corrections 
but it is argued that the qualitative behavior should be 
reliable. We merely note here that the dependence of 
eq. ([T]) on the constituent quark mass is rather prob- 
lematic. In the first place, the notion of a constituent 
quark mass is a concept which makes sense in the con- 
text of models but is not well defined within QCD. More- 
over, the concept becomes particularly problematic in the 
strong field regime where the typical lengths probed are 
short — shorter than the characteristic size expected of a 
constituent quark. 



Given this situation, it is important to see whether 
there is a viable model-independent way to study the 
strong magnetic field case. We argue in this paper 
that one can apply the techniques of perturbative QCD 
(pQCD) straightforwardly to this system, provided that 
one accounts for the strong fields by including the effect 
of the field on quark propagators to all orders in the field 
strength — which can be done using a proper time formal- 
ism. The central insight here is that when qB ^ A^ the 
characteristic Landau orbit size for a quark is small on 
the scale of A, and this sets the scale for the coupling 
when the quark interacts via gluon exchange with other 
quarks. Due to asymptotic freedom, the quark dynamics 
in the large field region corresponds to weakly coupled 
quarks; the principal way the vacuum responds to the 
field is simply to the rearrangement of states relatively 
deep in the Dirac sea due to the fields. In essence, when 
the characteristic size of a Landau orbit is much smaller 
than the typical size of a hadron, the quark in responding 
to the magnetic field doesn't "know" about the existence 
of hadrons and acts to very good approximation as a 
free particle. Corrections to this picture due to gluon 
exchange can be handled straightforwardly using pertur- 
bative methods. 

We focus on the magnetization as our probe of vac- 
uum response. In the next section we compute the mag- 
netization at leading order (one-loop) perturbative cal- 
culation. Following this we show that at next-to-leading 
order (two-loops) there is no correction to the one-loop 
result. We end with a brief discussion of the implications 
of these calculations. 



II. MAGNETIZATION AT LEADING ORDER 

The magnetization can be calculated as the derivative 
of the vacuum energy with respect to the magnetic field. 
In general, the magnetization can be calculated as [13] 
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calculated using the lagrangian appropriate 
to our assumptions. 

Here we consider the case of eB Aqc d so the nat- 
ural way to compute the effective Lagrangian pQCD re- 
quires the relevant scales in the problem to be much 
larger than Kqcd- 

The lowest-order effective Lagrangian is obtained from 
QCD at the one-quark loop level. The logic here is simply 
that in pQCD we are doing an expansion in Ois which we 
take to be small at these scales. At this order, the various 
flavors of quark do not communicate with each other — 
the up quark contribution to the effective exaction at the 
one-loop level is completely independent of down quark 
properties (charge and mass). It is only at three-loop 
order or higher that one has contributions from more 
than a single flavor. Thus at this order one can simply 
compute the contribution from each flavor independently. 
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We note here that the analysis depends on the mass of 
the quarks being smaU. However, in this context smaU 
merely means that ^ eB. Thus for sufficiently large 
fields the formalism applies even to heavy quarks. In 
intermediate regions, say, eB ~ 2GeV^, one can compute 
the perturbative contribution to the magnetization from 
the u,d and s quarks; in this region contributions from 
heavier quarks will be suppressed. 

The contribution to the effective Lagrangian at this or- 
der from a given flavor of quark, /, is that of free fermions 
in the presence of a constant external magnetic field; thus 
it is formally identical to that computed in Schwinger's 
classic paper namely: 
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where superscript indicates the contribution to the 
effective action for a particular flavor of quark, qj is the 
charge for quarks of that flavor, and ruq is the current 
quark mass for that flavor. Changing variables and tak- 
ing the derivative to get the contribution to the magne- 
tization, we find that at 
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where the subscript (0) indicates that it is lowest order 
in as- The 2/3z term in the parenthesis has its ori- 
gin in the renormalization of the electromagnetic field 
to due fermion loops; its presence ensures that the inte- 
grand converges as z ^ 0. Because B ^ m^, the expo- 
nential decays very slowly in the integrand of I{B,mq); 
thus, the large z region dominates. It is then conve- 
nient to separate out the large z contribution; we do this 
by dividing the integral into high z and low z contribu- 
tions divided by some number c with the properties that 

O l,^c< 1: 
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In the region of interest, I2 dominates. Moreover, in this 
region one can drop the coth(z) and z/ sinh^(z) terms in 
J2, since in this region the -2/3z: dominates. Thus, up to 
small corrections. 



IiB,m,) = -\j^'^-^e-</^^'^^^. 
Evaluating the integral we find that I{B, niq 
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r(°) [Z) is a plica function of zero order. In the region of 
interest ( — g-c ^ 1), this reduces to: 
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up to small corrections, where /i^ is a scale parameter. 
The leftover constant in eq. ^ is of the same order as 
previously neglected terms and can be dropped. Putting 
this together yields 
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up to small corrections. The leading corrections to this 
are proportional to B itself and, by a judicious choice of 
/i, these can be fully canceled. Since important contri- 
butions to these canceled terms come from the nonpcr- 
turbative region, /i^ encodes important nonperturbative 
physics. Note, however, that in the limit of very strong 
fields, the final answer is very insensitive to the precise 
value of fi. 

Equation ([S]) is a principal result of this paper. 



III. MAGNETIZATION AT NEXT-TO-LEADING 
ORDER 

Having determined the leading order perturbative con- 
tribution to the magnetization, it is important to consider 
corrections. These may be of two sorts — perturbative 
and nonperturbative. The nonperturbative corrections 
come from low momentum physics and as such cannot be 
computed directly from QCD via presently known ana- 
lytic techniques. However, such effect should be power- 
law suppressed in Aqq^^/ (eB) compared to the leading- 
order result of eq. ([5]) and thus very small at large fields. 
One expects the dominant corrections to eq. ([8]) to be 
perturbative. Such corrections to the magnetization can 
be written in the form 



(9) 



where the c's are dimensionless constants and as, the 
strong coupling constant, is evaluated at the scale of the 
problem: {eB). In this section, we evaluate the leading 
pertubative correction — ci. We show that it vanishes. 

The calculation proceeds from the effective action at 
next-to-leading order. This corresponds to the vacuum 
diagram depicted in fig. [TJ Note, as with the leading- 
order effective actions, the contributions from each fiavor 
are isolated from each other — each diagram contains only 
one flavor of quark. Thus at this order one can still com- 
pute the contribution from each flavor separately. The 
rules for calculating this diagram are identical to the 
usual QCD Feynman rules with one important difference: 
the quark propagator is replaced by the propagator for a 
spin-i fermion in a constant magnetic fieldflG], and then 
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the expression for the contribution to the effective La- 
grangian at next-to-leading order for a given flavor can 
be written as follows: 



FIG. 1: Diagram contributing to magnetization at next-to- 
leading order 
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where the superscript "1" indicates an expression at or- 
der ttg. 

Note that the expression for the propogator c;ontaiiis 
the quark mass. In the high field limit, perturbative ex- 
pressions should not depend on the quark mass. We in- 
clude it here merely to serve as an infrared regulator. 

The integral I{s,s',B,mq) is divergent as written. 
However, we are only interested in the B dependence. 
We can subtract off a divergent constant — the value at 
B — without affecting the magnetization. The in- 
dependent part of the total then becomes 

POO POO 
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This integral is convergent; however, it cannot be ex- 
pressed analytically in closed form. Fortunately, we are 
interested in the regime of eB ^ rrig. Thus, we can im- 
plement the same trick as was employed for the one-loop 
expression: divide the integral into a high and low s re- 
gion with the knowledge that the dominant contribution 



comes from a large s region: we obtain 
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up to small corrections, where mg is the separation scale 
introduced. Working in the regime juq 3> mg and taking 
the derivative with B, we find a dependence of the form 
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where is some renormalization scale chosen to minimize 
corrections due to nonperturbative effects. Combining 
this with the lowest-order result, yields 
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plus corrections of order a^. 

The critical thing to notice is that the correction term 
in eq. (|14p is proportional to asiri^/ (qfB) times the 
lowest-order term rather than just a constant times as- 
Thus, it is power-law suppressed in B and the coefhcient 
ci is zero. In fact, for light quarks this correction term 
is simply not reliable — it represents a small piece of the 
nonperturbative power-suppressed contribution and can- 
not be separated from it in a meaningful way. For the 
case of heavy quarks with very strong fields (the regime 
qjB rrig ^qcd) correction term is presumably 
reliable: it still represents a nonperturbative correction 
that is power-law suppressed at strong fields, but it is 
the leading contribution to this in a heavy quark expan- 
sion. Regardless of whether the quark masses are small 
or large, eq. ^ is accurate up to corrections of order aj.. 



IV. CONCLUSIONS 

We have examined the magnetization of the QCD vac- 
uum in the presence of a constant magnetic field in the 
strong field regime, eB ^ Aqcd, where perturbative 
QCD should be valid. We have found the contribution 
from any fiavor of quark with qfB 3> rri^. Combining 
these together, one finds that, in a regime in which the 
eB is either much greater or much smaller than all of the 
quark masses, the magnetization can be written as 
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where active flavors refer to those whose squared mass is 
well below eB. 

One striking feature of this result is that the two- 
loop correction makes no perturbative contribution. At 
present, we are unsure as to whether this indicates some- 
thing deep about the underlying structure of the theory, 
or rather it is something of an accident. 

Much of the motivation for this work was to obtain a 
model-independent prediction for the strong field region. 
Thus it is important to explore the relationship of our 
result to eq. ^ obtained in ref. [l3 using a plausible 



model and heuristic reasoning. At first glance, the two 
results are rather similar once one recognizes that the 
constituent quark mass is of 0(A) in eq. ([1]): both calcu- 
lations have the magnetization growing linearly with B 
times a slowly increasing function of B/A — in the case of 
eq. p5|) . logarithmically (noting that /i A), and in the 
case of eq. ([!]), a very small power law. Thus, it might 
appear that these are two essentially complementary de- 
scriptions of essentially the same physics. However, when 
one looks closely, there are important qualitative differ- 
ences between the two descriptions. 

From a theoretical perspective, one can use the number 
of colors as a probe of the nature of the dvnamics fiol. [20l| . 
Note that the perturbative expression derived in eq. (fT^ 
has a magnetization proportional to Nc- The origin of 
this is clear. The active players in dynamics are quarks — 
which act nearly as free particles for the purposes of com- 
puting the magnetization, since the length scale of Lan- 
dau orbits is small compared to characteristic hadronic 
sizes. Since the number of quarks participating is propor- 
tional to Nc, so is the magnetization. In contrast, eq. ([1]) 
has no overall factor of Nc- This reflects the dynamics 
modeled in ref. [12], where the magnetization was driven 
by a condensation of color-singlet "mesons" . Moreover, 
phenomenologically, the two results differ qualitatively 
in the regime of extremely large fields. In such a regime 
the power-law growth in eq. ([T]) ultimately becomes much 
larger than the logarithm in eq. (fT5|l . 

We believe that in the extreme strong field regime, 
the perturbative calculation described here should be- 
come increasingly accurate for reasons discussed above: 
at such extreme fields the localization scale of a Landau 
orbit is far smaller than the size of hadrons and, as far 
as vacuum polarization effects are concerned, should ap- 
pear to be essentially free. Thus, we suspect that the 
model described in ref. [12] must break down by the time 
the fields get very strong. This is highly plausible. The 
model depends on the notion of constituent quarks. This 
concept, while of possible utility in describing low- lying 
hadrons, is not appropriate for describing the perturba- 
tive regime. 

Support of the U. S. Department of Energy under 
grant number DE-FG02-93ER-40762 is gratefully ac- 
knowledged. 



[1] For a review, see: C. F. Perdrisat, V. Punjabi, and 

M. Vanderhaeghen, Prog. Part. Nucl. Phys.59 (2007)694- 

764, hep-ph/0612014 
[2] For a review, see: M. Vanderhaeghen, Eur. Phys. J. A8 

(2000) 455-479. 

[3] ELI project, http:/ /www.extreme-light-infrastructure.eu/ 
[4] D. AUor, P. Bedaque, T. D. Cohen, Phys. Rev. C 75 

(2007) 034001 

[5] A. Yu. Babansky, E. V. Gorbar, G. V. Shchepanyuk, 

Phys. Lett. B419 (1998) 272 
[6] S. P. Klevansky and R. H. Lemmer, Phys. Rev. D 39 



(1989) 3478 

[7] K. G. Klimenko, Z. Phys. C54 (1992) 323; K. G. Kli- 
menko, hep-ph/9809218 ; D. Ebert, K. G. Klimenko, 
M. A. Vdovichenko, and A. S. Vshivtsev, Phys. Rev. D 
61 (1999) 025005, hep-ph/9905253 

[8] V. P. Gusynin, V. A. Miransky, I. A. Shovkovy, 
Phys. Rev. Lett. 73 (1994) 3499; V. P. Gusynin, 
V. A. Miransky, I. A. Shovkovy, Phys. Lett. B349 (1995) 
477 

[9] A. Goyal and M. Dahiya, Phys. Rev. D 62 (2000) 025022 
[10] H. Sugamura and T. Tatsumi, Ann. Phys. 208 (1991) 



6 



470; 

[11] S. Schramm, B. MuUer, and A. J. Schramm, 
Mod. Phys. Lett. A7 (1992) 973; 

[12] D. Kabat, K. Lee, and E. Weinberg, 

[13] T. D. Cohen and R. F. Lebed, Phys. Rev. D 74 (2006) 
056006, hep-ph/0608038 ; T. D. Cohen and W. Bro- 
niowski, Phys. Lett. B292 (1992) 5, hep-ph/9208253 
T. D. Cohen, Phys. Lett. B359 (1995) 23, 
hep-ph/9505427; T. D. Cohen, Phys. Lett. B554 (2003) 
28, hep-ph/0210278; T. D. Cohen, Rev. Mod. Phys. 68 

(1996) 599 

[14] L A. Shushpanov and A. V. Smilga, Phys. Lett. B402 

(1997) 351 



[15] J. Gasser and H. Leutwyler, Annals Phys. 158 
(1984) 142; J. Gasser and H. Leutwyler, Nucl. Phys. 
B250 (1985) 465; for reviews, see G. Ecker, 
Prog. Part. Nucl. Phys. 35 (1995) 1; A. Pich, 
Rept. Prog. Phys. 58 (1995) 563 

[16] J. Schwinger, Phys. Rev. 82 (1951) 664 

[17] T. D. Cohen, D. A. McGady and E. S. Werbos, 
Phys. Rev. C 76 (2007) 055201 

[18] E. S. Werbos, Phys. Rev. C 77 (2008) 065202 

[19] E. Witten, NucL Phys. B160 (1979) 57 

[20] G. 't Hooft, Nucl. Phys. B72 (1974) 461; G. 't Hooft, 
Nucl. Phys. B75 (1974) 461 



